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ABSTRACT 
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When the 2dF QSO Redshift Survey (2QZ) is complete, a powerful geometric 
■ test for the cosmological constant will be available. By comparing the clustering along 

\D 1 and across the line of sight and modelling the effects of peculiar velocities and bulk 

motions in redshift space, geometric distortions, which occur if the wrong cosmology 
is assumed, can be detected. 

In this paper we investigate the effect of geometric and redshift-space distor- 
<i ■ tions in the power spectrum parallel and perpendicular to the observer's line of sight, 

Qh P S (k\\ , kj_). Ballinger et al. developed a model to estimate the cosmological constant, 

q A, and the important parameter (3 « O^ 6 /b from these distortions. We apply this 

$_i ' model to a detailed simulation of the final 25k 2QZ, produced using the Virgo Con- 

\ sortium's huge Hubble Volume N-body A-CDM light cone simulation. We confirm the 

& ■ conclusions of Ballinger et al. ; the shape of the redshift-space and geometric distortions 

are very similar, and discriminating between the two to produce a purely geometric 
constraint on A is difficult. When all the uncertainties in measuring P s (fc||,kjJ for 
the 2QZ are taken into account we find that only a joint A — (3 constraint is possible. 

By combining this result with a second constraint based on mass clustering evolu- 
tion, however, we can make significant progress. We predict that this method should 
allow us to constrain (3 to approximately ±0.1, and fl m to ±0.25 using the final cat- 
alogue. We apply the method to the 2QZ catalogue of 10000 QSOs and find that 
this incomplete catalogue marginally favours a A cosmology, obtaining best fit val- 
ues of /3 = 0.39±g:if and Q. m = 1 - fi A = 0.231^. However, Einstein - de Sitter 
(f2 m = 1.0, fl\ = 0.0) models are only rejected at the lAa level in the current survey. 
The rejection of lambda-dominated (f2 m = 0.0, fl\ = 1.0) models is stronger at ~ 2a. 
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1 INTRODUCTION 



The large-scale structure of the Universe represents one of 
the most powerful discriminants between cosmological mod- 
els. We are using the AAT 2dF facility to make a redshift 
survey of some 25000 bj < 20.85™ QSOs (Boyle et al. 2000, 
Croom et al. 2001a,b). QSOs are highly effective probes of 
the structure of the Universe over a wide range of scales and 



can trace clustering evolution over a look-back time which 
is 70-80% of its present age. 

The 2dF QSO Redshift Survey (2QZ) comprises two 
5 x 75 deg 2 declination strips, one at the South Galactic 
Pole and one in an equatorial region in the North Galac- 
tic Cap. QSOs are selected by ultra-violet excess (UVX) 
in the U-B-.B-R plane, using APM measurements of UK 
Schmidt Telescope photographic plates (Smith et al. 2001). 
The 2dF spectra obtained are identified automatically by 
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AUTOZ (Miller et al., in preparation), which also determines 
the redshift of the QSOs. These identifications are confirmed 
by two independent observers. Approximately 18000 QSO 
redshifts have been obtained to date, making this survey al- 
ready an order of magnitude larger than any previous QSO 
survey. 

Observations of high redshift Type la Supernovae 
(SNIa) have recently been used to measure the cosmological 
constant, A (e.g. Perlmutter et al. 1999; Reiss et al. 1998). 
By making the simple assumption that SNIa are standard 
candles, one can use the apparent luminosity and redshift of 
such supernovae to constrain the distance-redshift relation, 
hence determining the value of A. We currently don't have 
a good understanding of SNIa physical processes, and there 
are still some uncertainties regarding the standard candle 
assumption. For example, the above analyses apply a lit- 
tle understood empirical correction to the peak magnitude 
based on the decline rate of the SN light curve. There have 
also been suggestions that the SN light curves may be af- 
fected by cosmic dust (e.g. Aguirre 1999), or evolve with 
redshift (e.g. Drell, Loredo & Wasserman 2000; Reiss et al. 
1999a,b; see however Aldering, Knop & Nugent 2000). 

Alcock & Paczynski (1979) suggested that A might also 
be measured from distortions in the shape of large-scale 
structure. Geometrical distortions occur if the wrong cos- 
mology is assumed, due to the different dependence on cos- 
mology of the redshift-distance relation along and across the 
line of sight. Phillipps (1994) suggested that by making the 
simple assumption that clustering in real-space is on aver- 
age spherically symmetric, geometrical distortions due to A 
could be detected by comparing the distribution of radial 
and transverse separations of QSO pairs. This approach has 
the advantage that it depends only on the assumption of 
isotropy, without the need for a standard candle. We apply 
a related test to the 2QZ QSOs. We can measure the cluster- 
ing of QSOs parallel and perpendicular to the observer's line 
of sight, through the two-point correlation function, or in 
this case, the power spectrum, P S (k^ , kj_). Only if we make 
the correct assumption about geometry would no anisotropic 
geometric distortion be detected. 

Unfortunately the problem is complicated by the fact 
that redshift-space distortions are also caused by peculiar 
velocities. On the large, linear scales probed by the power 
spectrum coherent peculiar velocities due to infall into over- 
dense regions are the main cause of anisotropy. According 
to gravitational instability theory these distortions depend 
only on the density and bias parameters via /3 ~ Q^ 6 /b. 

In this paper we investigate the redshift-space distor- 
tions in the power spectrum parallel and perpendicular to 
the observer's line of sight, P s (k^,k.±), using mock 2QZ 
catalogues generated from the Hubble Volume simulation 
(Frenk et al. 2000). The aim is to determine how well we 
will be able to measure the important cosmological param- 
eters A and (5 using the 2QZ observations. 

We apply the model from Ballinger et al. (1996) for the 
redshift-space distortions to a measurement of P s (fc||,kj_) 
from the simulated QSO catalogue, and derive a joint con- 
straint on A and (5. The angular dependence of geometrical 
distortions differs only slightly from that of the redshift- 
space distortions, and so the constraint is fairly degenerate. 
To break the degeneracy, different constraints on A and j3 
need to be considered. Here we examine one possibility based 



on the evolution of mass clustering from the average redshift 
of the 2QZ to the present day. Finally we apply the method 
to the first public release of the 2QZ; the 10k Catalogue 
(Croom et al. 2001b). 



2 POWER SPECTRUM ANALYSIS 
2.1 Mock 2dF QSO Catalogues 

In order to test how effectively the 2QZ will measure large- 
scale structure, we are using the Virgo Consortium's huge 
Hubble Volume N-body ACDM simulation (Frenk et al. 
2000) to simulate three light cone strips of the 2QZ. The 
parameters of the simulation are f2b=0.04, f2cDM=0.26, 
S7a=0.7, H o —70 km s _1 Mpc _1 and the normalisation, as, 
is 0.9, consistent with the abundance of hot X-ray clus- 
ters (White, Efstathiou & Frenk 1993) and with the level 
of anisotropies in the cosmic microwave background radi- 
ation (CMB) found by COBE (Smoot et al. 1992). One 
billion mass particles are contained within a cube that is 
3,000/i~ Mpc on a side. One of the vertices was chosen to 
be the observer and the long axis of the light-cone was ori- 
ented along the maximal diagonal. The light-cone therefore 
extends to a depth of ~5,000/i _1 Mpc, which corresponds to 
z ~4 in the ACDM cosmology. The solid angle of the light- 
cone is 75x15 degrees which is split into three 75x5 degree 
slices. The 2QZ consists of 2 such slices, and so to generate 
the mock catalogue we use the two disconnected slices. We 
also consider the third slice for the error analysis. 

To create realistic mock catalogues, we must bias the 
mass particles to give a similar clustering pattern, and the 
same angular and radial selection function and sampling 
expected in the final 2dF Survey. The biasing prescription 
adopted follows method 2 as described by Cole et al. (1998), 
with the parameters chosen to produce a constant clustering 
amplitude with redshift - as seen in previous surveys (Croom 
& Shanks 1996). It is an Eulerian, or final density scheme 
due to the nature of the N-body simulation output, and is 
based only on the local density field, implying asymptotic 
scale-independence of the bias on the large scales we are 
probing. The biased particles were then sparsely sampled to 
match the radial distribution of the QSOs in the preliminary 
2dF catalogue, and the expected number density of the final 
catalogue over the redshift range 0.3 < z < 2.2. For further 
details see Hoyle (2000). 

The Hubble Volume simulation was designed to probe 
large-scale structure, at the expense of spatial resolution. 
As a result the simulation does not probe any structure on 
scales 3 h" 1 Mpc. In this analysis we are binning the data 
onto a grid with bin size ~ 12 h~ x Mpc, and only considering 
power on scales 32 ft -1 Mpc, so we expect this to have 
little effect on our final results. 



2.2 Power spectrum estimation 

The power spectrum estimation is carried out as described in 
Outram, Hoyle & Shanks (2001), using the method outlined 
in Hoyle et al. (1999) and Tadros & Efstathiou (1996). In 
orde r to apply the distant observer approximation (see sec- 
tion 2.3) the mock QSO survey is divided into four regions. 
Each region is embedded into a larger cubical volume, which 
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Figure 1. P s (fcy , kx) estimated from the 25k mock QSO survey. 
Filled contours of constant log(P(fc)/h _3 Mpc 3 ) are shown as a 
function of fcy/hMpc - 1 and kx/^Mpc - . Overlaid is the best- 
fit model (see section for P s (k^ , kx), with a fiat C m = 0.33 
universe, with f) = 0.33 and a v = 400 km s . 

is rotated such that the central line of sight lies along the 
same axis of the cube in each case. The density field is binned 
onto a 256 3 mesh, using nearest grid-point assignment. The 
power spectrum of each region is estimated using a Fast 
Fourier Transform (FFT), and the average of the resulting 
power spectra is taken. The result, binned logarithmically 
into A: 1 1 and kx, is plotted in figure jjj. 

2.3 Survey geometry 

The 2QZ covers two declination strips of approximately 
5° x 75°, and we have three such mock declination strips 
available. When calculating P s (fc||,kx) information about 
the line of sight must be retained. To achieve this the data 
are divided into subsamples that subtend a small solid angle 
on the sky, and the distant-observer approximation is then 
applied. The data from each declination strip are split into 
4 regions each approximately 5° x 20° during this analysis, 
and the average of the resulting power spectra is taken. 

This approach introduces a small systematic error; the 
measured P s (fc||,kx) is in fact a convolution between the 
true power spectrum and a function related to the sample 
window introduced due to the assumption of parallel lines of 
sight. The resulting power spectrum is more isotropic, and 
this leads to a small, yet systematic under-estimate of /3. 
For a 5° x 20° window this error is < 1% and considerably 
smaller than the statistical errors obtained. Although it is 
possible to correct for this effect (see Cole et al. 1994), for 
the purposes of this analysis we have chosen to ignore it. 

The measured power spectrum is also convolved with 



Distortions in the 2dF QSO Power Spectrum 3 




0.01 0.10 

Figure 2. The power spectrum of the 25k mock QSO survey 
window function. Filled contours of constant log(| W\ 2 /h~' i Mpc 3 ) 
are shown as a function of fcy/hMpc -1 and kx//iMpc -1 . To 
minimise the effects of the window function, we consider only 
the region bound by the white line above, with fcn & |kx| > 
0.02 h Mpc^ 1 . 



the power spectrum of the window function. As the QSO 
co-moving number density is almost constant out to very 
large scales (Hoyle et al. 2001a), we are effectively consider- 
ing a volume limited sample. Hence each QSO carries equal 
weight, and the survey window function, W(x), simply takes 
a value of unity in the volume of the universe included in 
the survey, and zero elsewhere. This is approximated using 
a catalogue containing a large number of unclustered points 
with the same radial and angular distribution as the survey, 
and its power spectrum is calculated in a similar manner to 
that of the data. Figure ^ shows the power spectrum of the 
window function, \W\ 2 , as a function of fen and kx, and fig- 
ure ^ shows the window function plotted as a function of k, 
with a comparison against that of a normal P(k) analysis, 
binning in shells of k. 

The window function is sharply peaked in k-space, vary- 
ing as a steep power law, ~ k~ 3 or steeper; much steeper 
than the expected QSO power spectrum. Therefore the con- 
volution should have little effect at wavenumbers where W\ 2 
takes a value much less than about a tenth of its maximum 
value. There is some evidence in figure |^ that the window 
function is slightly shallower along, and in particular across 
the line of sight. The directional dependence of the window 
function also explains the scatter seen in figure This ex- 
tension at low values of fey is due to the relatively narrow 
width of the declination strips; 5° at z = 2.4 corresponds 
to r ~ 300ft _1 Mpc, or k ~ 0.02. To avoid any problems 
caused by such effects we consider only those regions with 
fc|, & |k ± | > 0.02 h Mpc -1 . 
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Figure 3. The power spectrum of the 25k mock QSO survey 
window function (points), shown as a function of k. The bold 
points are those which were used in the analysis (see Fig. [2]) , with 
wavenumbers with feii > 0.02, |k ± | > 0.02, and k < 0.2. The win- 
dow function is a steep power law, varying as ~ k~ 3 . The overlaid 
line is the window function from the usual power spectrum anal- 
ysis, binning in shells of fc, of the 25k mock catalogue. 



2.4 Error analysis 

The errors in the power spectrum are estimated in two ways. 
Firstly we have three separate strips on which the analy- 
sis is carried out, and we can estimate the error from the 
dispersion in the three measurements. Alternatively we can 
consider the error estimate of Feldman, Kaiser & Peacock 
(1994; FKP). A comparison of the two error estimates is 
shown in figure ^. The two estimates agree fairly well, but 
with large scatter. Although some of the scatter in the error 
estimates is intrinsic, the dependence of the error on the an- 
gle of the wavevector to the line of sight, as seen in figure || 
also contributes to the scatter. The FKP error is, on average, 
slightly lower than the error estimated from the dispersion 
across the mock strips. Neither approach is ideal; the error 
estimated from the dispersion has a high uncertainty due 
to the small number of mock catalogues, whereas the FKP 
error estimate assumes a thin shell geometry, and can under- 
estimate the error on some scales, due, for example, to the 
onset of non-linearity. To minimise such problems, we took 
the error in each bin to be the larger of the two estimates. 
The adopted fractional errors are shown in figure ^. 



3 MODELLING REDSHIFT DISTORTIONS 

The power spectrum model incorporating redshift distor- 
tions that we fit to the data is presented in Ballinger, Pea- 



Figure 4. A comparison of the fractional errors using FKP esti- 
mation (circles) and the dispersion across 3 mock strips (crosses). 
The dotted and solid lines show the average error, binned in shells 
of k, calculated using FKP estimation and the dispersion across 
3 mock strips respectively. 



- 0.60 . ^^fc 




0.1 

k L 

Figure 5. The adopted fractional errors in P s (k^,U.±) (the 
larger of the FKP error and the dispersion across the three cat- 
alogues in each bin) are shown as a function of fcy/ZiMpc -1 and 
kj^/hMpc — 1 . The fractional errors vary from < 0.2 in the darkest 
region to ~ 0.8 in the lightest region. 
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cock and Heavens (1996). For clarity and convenience we 
review the main details below. 



3.1 Geometrical distortions 

The power spectrum analysis was carried out assuming an 
fi m =0.3, Qa=0.7 universe, matching that of the Hubble Vol- 
ume simulation. If the true geometry differs from this then 
our distance calculations will be wrong be a factor f± per- 
pendicular to the line of sight, and fu along the line of sight 
(as defined in Ballinger et al. 1996). Thus we can define a 
geometric flattening factor: 



F(z) 



A 

fx 



(1) 



(Ballinger et al. 1996). For an Einstein - de Sitter (fi m = 1.0, 
Qa = 0.0), and lambda-dominated (fl m = 0.0, Qa = 1.0) 
universe, F w 0.81 and 1.53 respectively at the mean depth 
of the survey [z « 1.4). The effect this has on the power 
spectrum is given by 



anisotrop 



ic(fc||,k x ) = 



-ftrue (A?) 

f ± +n F 



(2) 



(Ballinger et al. 1996) where n — kn/k, and n is the spectral 
index of the power spectrum. 



3.2 Redshift distortions 

On large, linear scales the main cause of anisotropy are co- 
herent peculiar velocities due to the infall of galaxies into 
overdense regions. This anisotropy takes a very simple form 
in redshift-space, depending only on the density and bias 



parameters via the combination (3 fts Q^ /b: 



P s (fc|, 



P R (k) + 



21 2 



(3) 



(Kaiser 1987) where P s , and P R refer to the redshift-space 
and real-space power spectra respectively. The Kaiser for- 
mula also assumes the distant-observer approximation, and 
hence is valid in this case due to the method of estimation 

Of P S (fc||,k X ) (§0. 

We expect the parameters b, Q m , and hence j3 to vary 
as a function of redshift. To keep the number of free param- 
eters to a minimum, however, we have chosen to fit a sin- 
gle average value of (3qso{z ~ 1.4) over the redshift range, 
0.3 < z < 2.2, probed by the QSOs. Large variations of (3 
over this range could introduce a systematic error into our 
results, however, the value of f3(z ~ 1.4) determined from 
the first 10000 2dF QSOs suggests th at an y variation with 
redshift is relatively small (see section 5.2). We hope to in- 
vestigate the evolution of (3 in more detail using the final 
2QZ sample of 25000 QSOs. 

Anisotropy in the power-spectrum is dominated by the 
galaxy velocity dispersions in virialized clusters on small, 
non-linear scales. This is modelled by introducing a damping 
term. The line of sight pairwise velocity (o- p )Q distribution is 
modelled using an exponential, which produces a better fit 
to both observations and N-body simulations than either a 



* In power spectra a v is implicitly divided by Ho and quoted in 
units /i -1 Mpc. (Ho = 100/t km s _1 Mpc -1 ) 



Gaussian or power law distribution (Davis & Peebles 1983; 
Fisher et al. 1994). Observational redshift determination er- 
rors may add a second, Gaussian component to the obs erved 
velocity dispersion of the 2QZ QSOs (see section 



5.1 



). An 

exponential is given by a Lorentzian factor in redshift-space: 



D [kfj,<jp] 



1 + 1 (kfj,a p ) 



(4) 



This model is only a good approximation in mildly non- 
linear situations. Cole et al. (1995) used numerical sim- 
ulations to show that it is effective for wavelengths A > 
15h~ Mpc, or k ^ 0.4; considerably smaller than the scales 
we are considering. At these wavenumbers the difference be- 
tween models assuming a Gaussian or exponential velocity 
dispersion is negligible (Ballinger et al. 1996). 

Combining these effects leads to the final model: 



P s (fe||,ka 



P R (k) 
f ± +n F 

1 + JU 



1 + /X 



13+1 
F 2 



F 2 



(•») 



(Ballinger et al. 1996) where a' p — o~ p /fn. 

There are several free parameters in the model; F, f3, a p , 
and P R (k), the underlying real-space power spectrum. To 
reduce the uncertainties on the parameters of interest, (3, &iF 
we make simple assumption about the other parameters. a v 
has relatively little effect on the power spectrum at large 
scales and so is not well constrained in this analysis. We 
have chosen to fix it at 400 km s _1 ; the measured value of 
the Hubble Volume simulation at z ~ 1.4 (Hoyle 2000). This 
is similar to recent measurements from local galaxy surveys 
(Outram et al. 2001, Peacock et al. 2001). Whilst it may 
be expected that a v falls with increasing redshift, the finite 
resolution of the Hubble Volume simulation, ~ 3 h~ l Mpc, 
introduces an additional velocity dispersion of similar or- 
der that keeps the observed value artificially high. A similar 
effect due to redshift determination errors affects the true 
QSO sample and is discussed later. To test any systemat- 
ics introduced by this assumption we repeat the analysis 
with o~p — 200 km s _1 and a p = 600 km s _1 . The effects 
can be seen in figure [| There is a small systematic shift in 
the best-fit (3 — S7 m relation that becomes insignificant at 
low values of Q m . A similar analysis of redshift-space dis- 
tortions considering £(<r, tt) (Hoyle et al. 2001b) is sensitive 
to smaller-scale clustering and so should constrain o v in the 
2QZ more effectively. 

We reconstruct a real-space power spectrum self- 
consistently from the measured mock QSO redshift-space 
power spectrum (Hoyle et al. 2001a), by inverting the above 
redshift distortion equations in each cosmology for an aver- 
age value of (i. Although the input real-space power spec- 
trum is available for the simulation, this is not the case for 
the 2QZ, and so we choose not to use it. Finally, we assume 
a flat cosmology, consistent with the recent CMB results 
from experiments such as Boomerang (de Bernardis et al. 
2000) and Maxima (Balbi et al. 2000). We choose to fit the 
variable Sl mj and fix Qa = 1 — ^m- 

Figure |7| illustrates how the fractional anisotropy in 
P s (fc||,kx) varies with cosmology. The ratio of two power 
spectrum models, with (3 = 0.39, fi m = 0.23, and (3 — 0.19, 
Qm = 1.0 (shown in Fig. H) are plotted. To discriminate be- 
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Figure 6. Filled contours of increasing \ 2 i n the O m — /3 plane for 
the 25k mock QSO survey, assuming u v = 400 km s" 1 . The solid 
line represents the best fit value of /3 for each £7 m and the l-a error 
is given by the dashed line. The dot-dot-dot-dash lines above and 
below the solid line show the best fit value of j3 for each f2 m with 
<jp = 200 km s _1 and a p = 600 km s _1 respectively. Overlaid 
is the best-fit (dot-dash) and \-a (dot) values of /3 determined 
using the method described in section | The joint best fit values 
obtained are /3 = 0.33 and O m = 0.33. The star marks the input 
Hubble Volume model values of j3 q = 0.35 and f2 m = 0.3. 



Figure 7. The ratio of two power spectrum models, /3 = 0.39, 
tt m = 0.23, and fi = 0.19, Q m = 1.0 (shown in Fig.||). This figure 
illustrates how the fractional anisotropy in P s (fey , k^) varies with 
cosmology. To constrain j3 and A, we need to detect variations of 
order two in the power spectrum anisotropy caused by geometrical 
and redshift distortions. 



is highly degenerate, and any values along the solid line are 
acceptable. The joint best fit value (see below) of f3 = 0.33 
and fi m = 0.33 has x 2 = 100.32, only marginally higher. 



tween these two models, and hence constrain j3 and A, we 
need to detect variations of order two in the power spectrum 
due to geometrical and redshift distortions. 

3.3 x 2 fitting 

We fit the model to the data using a \ 2 technique: 

2 \ ^ (^"observed Pmodel) /~\ 

X = ^ ~ 2 (6) 

°p(fc) 

The errors, <Jp(k), are the larger of the FKP error and the 
dispersion across the three catalogues in each bin. Only 
those wavenumbers with k« > 0.02, |kj_| > 0.02, and k < 0.2 
are used in the fit. The former constraints are applied to re- 
move the effects of the window function, and the latter be- 
cause the FFT is unreliable at smaller scales. The Nyquist 
frequency of the transform is k = 0.26. This constraint 
also prevents excessive non-linearity, where the model breaks 
down. 

The /3, and f2 m parameters in the model are adjusted 
until a minimum value of \ 2 IS obtained. Figure ^| shows 
X 2 contours in the f2 m - j3 plane. Nominally the best fit 
value is /3 = 0.02, and fi m = 0.05, with x 2 = 98.46 over 69 
degrees of freedom. This suggests that the errors have been 
slightly under-estimated. However, as figure shows, the fit 



4 MASS CLUSTERING EVOLUTION 

Applying the model from Ballinger et al. (1996) for the 
redshift-space distortions to a measurement of P s (fc||,kx) 
from the simulated QSO catalogue has given us a joint con- 
straint on f2 m and (3. The angular dependence of geometrical 
distortions differs only slightly from that of the redshift- 
space distortions, and so the constraint is fairly degenerate. 
To break this degeneracy, a different constraint on S7 m and 
f3 needs to be considered. Here we examine one possibility 
based on the evolution of mass clustering from the average 
redshift of the 2QZ to the present day. 

One of the early results from the 2dF Galaxy Redshift 
Survey is a measurment of (3 s (z ~ 0) = 0.43 ±0.07 (Peacock 
et al. 2001), using redshift-space distortions in the two-point 
galaxy correlation function. For each cosmology (again we 
only consider flat cosmologies) the value of the galaxy-mass 
bias can be found from j3 s , which in turn gives the value of 
the mass correlation function if the galaxy correlation func- 
tion at 2 = is known, as shown below. The evolution in the 
mass correlation function can be calculated for this cosmol- 
ogy, and hence by comparing the z ~ 1.4 to mass correlation 
function to that of the QSOs, an estimate of the QSO bias 
factor can be obtained. Finally, this can be used to derive 
an estimate of f3 q as a function of cosmology. One possible 
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caveat with this technique is that we are comparing values of 
f3 g and (3 q that were measured using different estimators and 
on different scales. To minimise any possible systematic ef- 
fects this analysis could be carried out self-consistently using 
the same method to determine /3 9 (Outram, Hoyle & Shanks 
2000). As the 2dF Galaxy Redshift Survey data is not yet 
available, we have chosen for simplicity to use a correlation 
function analysis. 

Rather than determining the value of the two-point cor- 
relation function at one particular point, we use the less 
noisy volume averaged two-point correlation function out to 
20/i -1 Mpc, £(20). To estimate the redshift space, volume 
averaged two-point correlation function of galaxies, we 
assume that the galaxy correlation function can be approx- 
imated by a power law of the form £| = (s/6/i. _1 Mpc)~ 1,7 . 
The power law approximation is in very good agreement 
with early results from the 2dF Galaxy Redshift Survey 
two-point correlation function over the range of scales 2< 
s <20/i _1 Mpc (Peder Norberg, private communication). 
Once the 2dF Galaxy Redshift Survey is completed, there 
will be no need to make this approximation. £ 9 is then found 
via 



€(20) = 



Jo 



4irs 2 ds 



= 0.3. 



(7) 



By integrating out to 20/i _1 Mpc, the non-linear effects on 
the volume averaged correlation function should be small, 
and are not considered in this analysis. Redshift measure- 
ment errors should only be a factor on small, ^ 5/i _1 Mpc, 
scales and so are also not considered. 

For each cosmology, the bias between the galaxies and 
the mass at z~0, b gp (0), can be found from 



V(o) 



»r(o) 

&(0) ' 



(8) 



The real space galaxy correlation function can be deter- 
mined from the redshift space galaxy correlation function 



f r = 



[l + |3 g (0) + i/3!(0)] : 



0) 



where the superscripts r and s indicate real and redshift 
space respectively. The real space mass correlation function 
at z=0 can now be found as 



WY 



(10) 



The real space mass correlation function evolves according 
to linear theory such that 



go) 



(ii) 



where G(z) is the growth factor, which depends on cos- 
mology, found from the formula of Carroll, Press & Turner 
(1992). When fi m =l, G(z) = (1 + z). 

We now relate the correlation function of the mass at 
z measured in real space to the amplitude of the QSO 
clustering at z, measured in redshift space as we wish to 
know /3 q (z), where the subscript q stands for QSO, as a 
function of Sim (0) . £|(20), measured over the redshift range 
0.3 < z < 2.2, is determined separately in each cosmological 
model from the latest 2QZ results (see Croom et al. (2001a) 



0.8 
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Figure 8. Significance contours given by joint consideration of 
the two constraints from the 25k Hubble Volume mock catalogues. 
The joint best fit values obtained are /3 = 0.33 and f2 m = 0.33. 



for details). In this analysis we have not taken into account 
the effect of redshift determination errors on the measure- 
ment of the QSO correlation function. This may lead to a 
slight systematic overestimate of the QSO bias. 
First we calculate fi m (z) using 

O m (0)(l + *) 3 

ilm[z) - n m (o)(i + z)3 + n A ( ' 

(valid for flat cosmologies only). P q (z) is then given by 
ft <z)°-° 

"'« = W' (13) 



where b qp (z) is defined by 



(,q(z) can be found via 



£(*) 



(14) 



(15) 



As £q is remeasured in each cosmological model, this is valid 
without consideration of geometric distortions. Combining 
equations [l3||l^ , and |l5| gives 



Qm(z 



i.2 m 



(16) 



a quadratic in fi'q(z) which can easily be solved, allowing 
/3 q (z) for any fi m (0) to be found. We substitute z = 1.4 to 
find the value of /3 q (z) at the average redshift of the survey. 

The errors on /3 q (z) are found in the standard way, i.e. 
by differentiating (3 q (z) with respect to all the variables that 
contribute to the error and summing the components of the 
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error in quadrature. In this case there are errors on /3 9 (0), 
£*(z) and £g(0). The error on (3 g (0) is ± 0.07, the error on 
£q(z) is the Poisson error found from the total DD counts. 
The error on £ 9 is the least well known as the power law 
approximation is used. We assume an error of around 20%, 
which is larger than the error expected from the completed 
2dF Galaxy Redshift Survey. 

To produce a fully self-consistent constraint on the 
mock catalogue data, we need to consider measurements of 
/3 9 (0) & Q for the galaxies and £, q for the QSOs from the 
same simulation, rather than real data. To do this properly 
would require producing a mock galaxy catalogue. Alterna- 
tively we can appeal to the fact that the value of = 1.4) 
in an f2 m =0.3, Qa=0.7 universe, obtained above, is in very 
good agreement with that of the Hubble Volume. Also, the 
Hubble Volume QSO mock catalogue was designed to have 
a correlation function that closely matches the observed 
2dF QSO clustering properties (Hoyle 2000). Hence the con- 
straint we derive below is still valid. 

The value of (3 q (z) and the one-sigma errors are plotted 
on figure []. Although totally degenerate with the value of 
fi m , this provides a different constraint to that from the 
redshift-space distortions. This was derived using signifi- 
cantly smaller scales than the P s (fc||,kx) analysis, and so 
we can treat the results as independent and hence combine 
the likelihoods, yielding a much stronger fit. The joint best 
fit values obtained are f3 q = 0.33±0.10 and fi m = 0.33±0.22. 
The one, two and three-sigma significance contours are plot- 
ted in figure ||. These results are consistent with the input 
values of j3 q — 0.35 and fi m = 0.3, and an Einstein-de Sitter 
model is ruled out in the Hubble Volume at over three-sigma. 




Figure 9. P s (fc|| , kj_) estimated from the 2QZ 10k Cat- 
alogue containing 10681 QSOs. Filled contours of constant 
\og(P(k)/h~ 3 Mpc 3 ) are shown as a function of fc||/hMpc — 1 and 
k^/ft Mpc -1 . Overlaid are the best fit model (solid contours) 
with P = 0.39 and Q m = 0.23, and a model with f2 m = 1.0 and 
/3 = 0.19 (dashed contours). 



5 THE 2QZ 10K CATALOGUE 

We now apply this analysis to the 2QZ 10k Catalogue 
(Croom et al. 2001b). This catalogue contains the most spec- 
troscopically complete 2dF fields observed prior to Novem- 
ber 2000. The colour selection technique used to create the 
input catalogue of QSO candidates for identification using 
the 2dF instrument gives high completeness (>90%) out to 
2 ~ 2.2 (Boyle et al. 2000). Adopting this as our high red- 
shift limit to ensure high completeness, we use a sample of 
8935 QSOs with redshift 0.3 < z < 2.2 in the following 
power spectral analysis. 

The measured QSO power spectrum is a convolution of 
the intrinsic pow er sp ectrum with that of the window func- 
tion (see section 2.3). As the survey is still incomplete, it 



currently has a patchy angular selection function and a com- 
pleteness map has to be constructed. This is used to create 
the random catalogue that approximates the complicated 
window function in the analysis. Extinction due to galac- 
tic dust was taken into account using the maps of Schlegel, 
Finkbeiner & Davis (1998). The construction of the window 
function is described in more detail in Hoyle et al. (2001a). 

There is a small anisotropy introduced by the minimum 
separation of fibres on the 2dF instrument. As the fibres can- 
not be placed closer together than 30", close pairs of QSOs 
cannot be observed simultaneously, and this introduces an 
angular selection function against such pairs in the 10k cata- 
logue. This is being addressed through separate observations 
of the QSOs missed for the final catalogue. The scales we are 



probing are considerably larger than those affected by this 
selection effect; at z ~ 1.4, k ^ 0.2 corresponds to 6 <; 40'. 
Hence we expect that this selection function should have 
little effect on our results. 

Due to the lower number of QSOs, and a more compli- 
cated window function, the power spectrum measurements 
were very noisy, and only coarse binning in k« and kx was 
possible. Figure ^ shows P S {k\\ , kx) estimated from the 2QZ 
10k catalogue. The power spectrum of the window function 
of the 10k catalogue are shown in figure jilj Due to the cur- 
rent incomplete coverage, the points lie slightly above those 
of the mock catalogue. The window function is still sharply 
peaked in k-space, however, and so by considering only those 
regions with ku & |kj_ | > 0.02 h Mpc - , any window func- 
tion effects should still be negligible. 



5.1 The effect of redshift determination errors 

There is an uncertainty in determining QSO redshifts from 
low S/N spectra of Sz ~ 0.0035 (Croom et al. 2001b), which, 
at the average redshift of the survey, assuming ACDM, cor- 
responds to an uncertainty in the line-of-sight distance of 
~ 5 /i -1 Mpc. This introduces an apparent velocity disper- 
sion of o p ~ 600 km s -1 . By adding this in quadrature to 
the intrinsic small-scale velocity dispersion in the QSO pop- 
ulation, we therefore expect an observed velocity dispersion 
of o p ~ 600 — 800 km s _1 (assuming an intrinsic veloc- 
ity dispersion of o p ~ — 500 km s _1 ). A second possible 
source of redshift determination error, due to the intrinsic 
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variation in line centroids between QSOs is harder to es- 
timate. Although the intrinsic pairwise velocity dispersion 
is perhaps best modelled with an exponential, the appar- 
ent velocity dispersion introduced by redshift determination 
errors is more likely to have a Gaussian distribution, and 
hence ideally we should adopt a model that is a convolution 
of both. At the scales we are considering, however, there is 
virtually no difference between the two models (Ballinger et 
al. f996), and so we have continued to adopt an exponential 
model. 

To test this we introduced a 5 h^ 1 Mpc Gaussian dis- 
persion to the line-of-sight distances of the mock QSOs from 
the Hubble Volume simulation and repeated the analysis. 
Still assuming a p = 400 km s _1 , the best-fitting value of 
f2 m was system atica lly higher for any given value of j3. As 
shown in section 3.2 an under-estimation of a p would lead to 



such a systematic shift. This was removed by increasing the 
value of tip in the model to take into account the apparent 
velocity dispersion introduced due to redshift determination 
error, assuming instead a v alue of a p — 700 km s _1 . 

As discussed in section 3.2, and shown in figure [| a v is 
not well constrained in this analysis, and so we have chosen 
to fix it at a p = 700 km s _1 for the 2QZ analysis. A comple- 
mentary approach, using the two-point correlation function, 
£(<7, 7r) is more sensitive to non-linear effects, such as the ve- 
locity dispersion, as it probes smaller scales, and so should 
provide a stronger constraint on a v . A careful analysis of the 
errors in redshift determination of 2QZ QSOs using AUTOZ 
would also constrain our assumed a v further. 



5.2 



X 2 fitting 



As before, we reconstructed a real-space power spectrum 
self-consistently from the measured 2QZ QSO redshift-space 
power spectrum (Hoyle et al. 200fa), by inverting the above 
redshift distortion equations for an average value of fi. The 
result of fitting the redshift-space distortions model to the 
observed power spectrum, P S (k\\ , kx), is shown in figure [flj 
A nominal best fit of Q m = 0.46, and j3 — 0.58 is obtained, 
with x 2 = 30.74 over 3f degrees of freedom. The fit is highly 
degenerate, however, and little constraint on either parame- 
ter alone can be reached. When we jointly consider the f2 m , 
/3 constraint from mass clustering evolution, as discussed in 
section |i|, a ACDM cosmology is marginally favoured with 
joint best fit values of j3 = 0.39±g;i? and Q m = 0.23^^3 be- 
ing obtained (see figure |f2|). The quoted errors include sta- 
tistical errors and the systematic error due to a ±1 00 km s _1 
uncertainty in a p . If the error in determining QSO redshifts 
were over-estimated by Croom et al. (2001 b) then we would 
find slightly lower values for A and (5. 

This result is in line with the recent determination of 
A from high reshift supernovae. Perlmutter et al. (1999) 
obtained an estimate of Sl m = 0.28_"q'o8, assuming a flat 
cosmology. Whilst statistically stronger than the geometric 
method, SNla physics is still poorly understood and it is 
possible that evolutionary effects could introduce large sys- 
tematic errors (Drell et al. 2000). The method described in 
this paper provides an independent determination of A that 
relies only on the assumption of isotropy, and relatively well 
understood redshift distortion models. 

The value of /3 q {z ~ 1.4) obtained is almost identical 
to that of local galaxies; e.g. /3 g (z ~ 0) = 0.43 ± 0.07 de- 



o 
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Figure 10. The power spectrum of the window function of the 
2QZ 10k catalogue, shown as a function of k. The crosses show 
the points which were used in the 10k analysis, with wavenumbers 
with k\\ > 0.02, |kjj > 0.02, and k < 0.2. The window function 
is a steep power law, varying as ~ k~ 3 . The points show the 
power spectrum of the window function from the analysis of the 
2QZ 25k mock catalogue. The 10k catalogue points are slightly 
higher than those of the full mock catalogue due to the current 
incomplete coverage of the survey. The overlaid line is the window 
function from the usual power spectrum analysis, binning in shells 
of k, of the 25k mock catalogue. 



termined using the 2dF Galaxy Redshift Survey (Peacock 
et al. 2001). If QSOs trace average galaxy environments, 
as suggested by deep imaging around QSOs (e.g. Croom & 
Shanks 1999), measurements of QSO clustering (Croom et 
al. 2001a, Hoyle et al. 2001a), and the observation that most, 
if not all, major galaxies appear to host a supermassive black 
hole (Magorrian 1998), then ft appears to vary very little as 
a function of redshift. With a larger QSO sample, we should 
be able to investigate the evolution of f3 by measuring it in 
narrower redshift bins. The main caveat would be whether 
we really are comparing the same population of objects at 
different redshifts. 



5.3 Comparison with Hubble Volume simulation 

Using the Hubble Volume simulation, we have demonstrated 
that a significant constraint on the A and /3 parameters will 
be available once the QSO survey is complete. The current 
survey, however, has fewer QSOs and a complicated window 
function. To investigate if this window function has any sys- 
tematic effect on the 10k results we can again turn to the 
Hubble Volume. 

We have repeated the power spectral analysis of the 
Hubble Volume simulation, using the current 2QZ win- 
dow function when constructing a mock catalogue of 10000 
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Figure 11. Filled contours of increasing \ 2 i n the n m - f3 plane 
for fits to the 2QZ 10k catalogue. The solid line represents the 
best fit value of /3 for each f2 m and the l-cr statistical error is 
given by the dashed line. The dot-dot-dot-dash lines above and 
below the solid line show the best fit value of j3 for each tt m with 
Up = 600 km s _1 and a v = 800 km s — 1 respectively. Overlaid is 
the best-fit (dot-dash) and l-cr (dot) values of /3 determined using 
the mass clustering evolution method described in section [l] The 
joint best fit values obtained are j3 = 0.39 and f2 m = 0.23. 

QSOs. The results are shown in figure [l^. The constraint 
on the A and (3 parameters obtained from the 10k Hub- 
ble Volume simulation is very similar to that obtained from 
the 10k 2QZ catalogue, although with slightly larger uncer- 
tainty. The joint best fit values obtained, taking into account 
the mass clustering evolution method described in section ^J, 
are f3 = 0.42 and fi m = 0.20. Comparing the results derived 
from the 10k and 25k mock catalogues, there is a shift in the 
best fitting parameters to a slightly higher (5 for a given fi m 
when the incomplete window function is considered, how- 
ever the two constraints are entirely consistent within the 
errors, indicating that possible systematic effects due to the 
10k window function are fairly small. If this shift is system- 
atic then taking it into account would slightly weaken the 
significance of our A constraint from the 10k 2QZ catalogue. 



6 SUMMARY AND CONCLUSIONS 

When the 2QZ is complete, a powerful geometric test for 
the cosmological constant (Alcock & Paczynski 1979) will 
be available. The distance between a pair of QSOs mea- 
sured along the line of sight, calculated from the relative 
redshifts, has a different dependence on the cosmological 
parameters than the distance measured in the angular di- 
rection. By comparing the clustering along and across the 
line of sight and modelling the effects of peculiar velocities 



Figure 12. Significance contours given by joint consideration 
of the two constraints from the 2QZ 10k catalogue. The errors 
from the power spectral analysis include statistical errors and the 
systematic error due to a ±100 km s _1 uncertainty in a v . The 
joint best fit values obtained are fi = 0.39 and H m = 0.23. 



and bulk motions in redshift space, geometric and redshift- 
space distortions can be detected. 

In this paper we have examined the potential of this 
test by considering the anisotropy in the 2QZ power spec- 
trum. Ballinger et al. (1996) developed a model to estimate 
A and /3 from the geometric and redshift-space distortions 
seen in the power spectra of redshift surveys. We have ap- 
plied this model to a detailed simulation of the 2QZ, pro- 
duced using the Virgo Consortium's huge Hubble Volume 
N-body A-CDM light cone simulation. The results confirm 
the conclusions of Ballinger et al.; the shape of the redshift- 
space and geometric distortions are very similar, and dis- 
criminating between the two to produce a purely geometric 
constraint on A is difficult. When all the uncertainties in 
measuring P s (fciukjj for the 2QZ are taken into account 
we have found that only a joint A — f3 constraint is possible. 

However, by combining this result with a second con- 
straint based on mass clustering evolution, we can make sig- 
nificant progress. We predict that this should allow us to 
constrain (5 to approximately ±0.1, and Q m to ±0.25 in the 
final 25k survey. Note that throughout this paper we have 
assumed a flat cosmology. 

One possible caveat with this technique is that we are 
comparing values of f3 g and (5 q that were potentially mea- 
sured using different estimators and on different scales. To 
minimise any possible systematic effects the value of f3 g 
could be measured using the same method (Outram, Hoyle 
& Shanks 2000). Although at a much lower redshift, there 
will also still be a small cosmological dependence in the re- 
sult which should be taken into account. 
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Figure 13. Filled contours of increasing \ 2 i n the n m - f3 plane 
for fits to the Hubble Volume mock 10k 2QZ catalogue. The solid 
lines represent the best fit value of fi for each f2 m and the 1-a 
statistical errors. Overlaid are the best-fit and 1-<t (dot) values 
determined for the 2QZ 10k catalogue (Fig. |nj), and the best- fit 
and 1-a (dashed) values for the 25k mock catalogue (Fig. ^). Tak- 
ing into account the mass clustering evolution method described 
in section ^, the joint best fit values obtained are /3 = 0.42 and 
fim = 0.20, compared to f3 = 0.33 and n m = 0.33 for the 25k 
mock survey. 



Applying this analysis to the 2QZ 10k Catalogue we 
find that a ACDM cosmology is favoured with joint best 



fit values of f3 = 0.39Io;i? and fi m = 0.23Io*3 being ob- 
tained. An uncertainty in determining QSO redshifts from 
low S/N spectra of Sz ~ 0.0035 introduces an apparent ve- 
locity dispersion of a p ~ 600 km s . The quoted errors in- 
clude a systematic error due to a ±100 km s _1 uncertainty 
in a v . A careful analysis of the errors in redshift determina- 
tion of 2QZ QSOs using AUTOZ would help constrain this; if 
the errors in determining QSO redshifts are currently over- 
estimated then we would find slightly lower values for A and 
13. A second possible source of redshift determination error, 
due to the intrinsic variation in line centroids between QSOs, 
is harder to estimate. 

A complementary approach, using the two-point corre- 
lation function, £(<r, 7r), to investigate redshift-space distor- 
tions in QSO clustering can also be considered (Matsubara 
& Suto 1996). This approach probes smaller scales, and so 
is more sensitive to non-linear effects, such as the velocity 
dispersion. A second paper will investigate the application 
of this analysis to the 2QZ (Hoyle et al. 2001a). 
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